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Optimum Wavelength Selection for Multi-Wavelength Absolute Phase 

Measuring Sensors 



p P Towfcrs, CE Towers, J D C Jones 

SchtmlofEfigutesrinffimdFfQfsicalScienDej, Ueriot-Wali Unfyetsf^, Edmbur^ £HI4 UK. 



1.0 Background 

Th.earc exists a multitade of Int^rf^pmetriic ba$ed ^mors 
for measming the phase of dectromagcietic wayes. 
Measutanib iSiat: effect the phase of a 'wsve include; 
defofmation [l], vibration [2], refractive iadex variiatiojx 
due to density changes and a Jbion of i»rdmilqTiea for single 
point and whole ^eld profilometiy &om synthetic 
j^rcure radar (SAR) [3] to j&jjige projectjion [4], to ap 
mterferBnce pattern ihere are two elements that contribiite 
TO die phase measurament dynamic xanga: the sub-&inge 
phase resohitlan and the number of fi'bises (fringe 
orders) spanned by tbe measuroaiesxt There is typically a 
simple Ainciion relating the measured optical fiinge 
phase to the desired measvrand. Tlxe process of phase 
measunxxg produces aub-i6TAge resolution, typically 
l/lOO'* to I/IOOO^** of a frmge. The sub-fringe resolution 
is calculated vi^ either pha^e steppxng [5] or Foarier 
transform techniques [6]. Howevsr, the interf gro meliic 
idiase is calculaled using an inverse trigonometo^ 
fiittctioD vrlth principal values over Ifae range -« to -Hr a± 
best Therefore the required phase Information is 
'wrapped' into an interval -with sharp discontinuities in 
the data at the edges of that interval which must be 
spatially or temporally unwrapped to obtain die ftinge 
order oofbimation. 

Z.0 Current Statc-of-the-Art 

A generi-c problem common to most mterferometnc 
techniques is the determination of absolute jfringe order 
in interttsrograms containing phase discontinuities or 
spatially or temporally discrete samples. Conventional 
interferometric analysis for single point data relies on a 
temporal scan to measure the t«lative pha&e diange in 
going from one state to another. In the case of whole 
field data, a spatial unwrapping is achieved using an 
appropriate algorithm giving relative information on the 
state (or ch9nge of state) across Qie image Said [7]. 

One well established research area for absolute phase 
ttteasutsemeixt is fiinge projection for profilometry wha;e 
the ^thetic wavelengths produced can be varied 

retetively easily. Therefore, techmqties ba^ed on 
tempoial phase unwrapping may be applied. Hutrtley et al 
ittttoduced a temporal phase unwrapping approach using 
a reverse exponential series of projected fringes given by: 
j,j-l,j-2,s-4,...»5/2 [8]. Here, the unwrappii:^ is 
perforraed between consecutive phase measurements to 
scale the fringe order calculated at the subsequent 



wavelength. Absolute fdnge order is obtained as the 
unwrapping is pcrformedl at each pixel independently 
along die time axis. Hie nmnber of projected ftiz]ge 
frsquenoies jreqiijred has been reduced to i$fy^^{s)'i-^. 
For each projected frhige freqv;ency fovir phase stepped 
images are obtained to determine the wrapped phase 
values in the interval -ic to -Hi. For exatnpl6, to measure 
over a range of 32 Binges 6 sets of 4 frames are needed, 
and &r 128 binges 8 set^ of 4 frames. Therefore, 
considerable thue is required to obtain the image frames, 
and ^ e^ch image may coqtiiin >3MB of irtfbtmatioin a 
aignificant data processmg problem is generated. 

3>0 Ba^is of tbe Inventiioni 

We have devised a novel sfxnftBgy for absolute fiiogo 
order identific&tioti in multi-wavelengdi heterodyne 
Interferometicy based csa optimum seledjon of due 
wavelengths to be used A theoredcal model of the 
process has been developed which allows the process 
reliability to be quantified. The methodology produces a 
Mwelength setcction which is optimum with respect to 
the minimum numbex of wavelengths required to achieve 
a target dynamic range. Cotkversely« the maximum 
dynamic muge is produced from a given number of 
optimcdiy selected wavelengths ntjiKsed in a sensor. The 
new concept introduced fbr optimum wavelength 
selection is scalable, i.e. from, a two wavelength system 
to a three waveletiglh system, iBom three wavelengths to 
foVDT etc. Therefore, whilst the previous technology 
allowed 2""* fringes to be counted absolutely where N is 
ihe number of fringe frequencies used, the new technique 
scales as wfaajre n is sin arfaitraxy real numfao: limited 
by phase noise of a practical interferometer. 

3*1 Theoretieal Development 

J././ Two Wavelength Inier^rometry 
hi two wavelength heterodyne inierferome<ry to 
eliminate $tep height ambiguity the difference in the 
number of fringes projected across the field of view must 
be £1. Let the number of projected fringes across the 
field be Nj^ fer the/* WJ^velengtix; then N^,^ -i^^i <1 

withiV^ii >//;ti« dijaerence in wrapped phases 

calculated at the two wavelengths, a heterodyne 
fbnctlon, can be expressed as a phase wiOnn the interval - 
n to 'it and will consist of a monotonic ramp across the 
Image. A convenient rej>resentatioD is to calculate a 
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discrete phase level for each of the fringe orders 
conesponding to iVju . This is obtained by gfubtractine a 
scaled version of the wrapped phase at //ji, &om ^bie 
heterodyne fiiaction - the scaling factor being given by 
i^Ai'^^izV^Ai- ^ practice, such a diacxete level 
heterodyne fionction can only identify a limited munber 
of fidnge ojdexa owhig to the presence of phai6 noi^e. 
JEadi Nji is known, as it is set in. a white light Bystegoi or 
can be measured in coherent systems [9]. 
Bach wrapped phase zneasurem^t confaixis phase noise, 
which is modelled bs a Gaussian distribution with zero 
meaxi and a standard deviation of cr^ [iO]. The 
hetcffodyne fiinctiotx then contaxo^ noise with a standard 
deviation given by •JZc^- We define a process 
robustness of 6<r, conespoadjrog to a probability of 
99.73% sncce^ in fringe order identification within &e 
mea^rement system. Therefore, the discrete p]tme teve]$ 

must be separated in phase by at least S^er^ . Hence, 
for 6ar reh^bility the numbBr of fringes which can be 
correctly identified is limited by; 

Defining the measurement dynamic range as the product 
of the pha^e resolution apd the number of fringes 
successfhlly numbered. Therefore, for a j^ase resolution 
of j/jOOthof a fringe^ 12 Imges can be nrnnbered to 
60- reliability, giving a dynamic range of 1200, which is 
insufficient for most engmeeixng applications. 

3 J, 2 Optimally Selected Wmelengths m iV 
Wavelength Inttirferomerty 

Ihe introduction of a third projected wavelength allows 
the generation of two independent heterodyne functions 
contaimng No^^ and iV^^ discrete levels* where 

•^wi^-^Dci =-^^1' discrete phase levels in each 
modified heterodyne jfunction wJU be seji^xiated byt 

, and S^^^ . (i) 

In this arrangement it is foimd that as increases 
i9^, decreases. Therefore, for nuodmum overall 
reliabiliiy in fringe order identification the optimimi set 
of projected fringe wavelengti]us is given fay the 
symmetrical arrangement where ^nn-^Ms an^A 
iViu, = Noli = V^ai • ^f wumbesT of discrete levels is 
not balanced the ejSfect is to increase one of Nj^j,^ and 
1 thereby bringing ^^^j or closer to the noise 
limit. Hence, from equation (l), fhe optimised S 

wavelength approach the c<mditlcin f^w: 6cr reliability 
fringe order identification is given b>% 



(2) 



The number of fringes which may be heterodyned 
reliably in an optimum 3 wavelength hctegrodyne setup is 
the;:efore the square of that for the two wavelengb case 
for the same phase measurement noise. 
Equations 1^ are used to define general expressions 
relating the numbers of projected frhiges at the N 
wavelengths. By assuming Nj^>Nf^ and 
' Nj^ , with 



W^Bii=^iW2"=—=^jD«ri optimum noise immunity, 
one general approach is giveaby, 

where: /i pumber of wavelengths, A^^ number of 
fringes in the largest fringe set, Nj, number of fringe^ in 
the i"^ fringe set Ihe tenn Ny, "0 is mcluded in order 
to generalize the expressions that follow. The general 
formolas to calculate the fringe order are givea by: 



(4) 
(5) 



fOT i=;l,....^;L-1 , where DL^ is tihe discrete level 

function, [DL^ is the r integer discrete level and H^j^ 

is the heterodyne between the 0 and wrapped phase 
map. The recursive jelationship for tlie integer discrete 

level funotions is hiitialtsed by setting IDL^ =0, The 
fringe order for the wr^ed phase map with Nj^ , fringes 
is given by IDL^ from equation (5). 
For example using four synthetic wavefengtihs with 
=64, then N^,^63, N^^^6Q, =48, fig, 1 
shows the discrete level functions resuWng from a 
simulation of the process. The phage dij6fe^ncE between 
discrete levels Cy^axis) is equal in each plot as expected 
for an opdnumi wavelength confrguration. A 
multiplicative intensity noise of has been applied to 
the cosine intensity frinses wJiich i? representative of 
scientifio CCD cameras fijod is a primary source of phase 
noise. 

Hie number of fringes that may be numbesred to 6<t 
reliability for a system with JL wavelengths is given by 
(from equation (2)) 



2jc 



The reliability of the process only depexids on the frrst 
and last -wrapped phase maps with all intcnnediate values 
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for other wavelengths cancelliog out (see equations 4 aad 
5). As the last wavelength corresponds to z»ro ^ing^^j 
jatfcer than perfonnlns a phase roeasufemeni here it is 
more accurate to set the wrapped pl^ase^ for ihia 
wavelength to zero. Hence the jifocess reliabiliiy depwds 
only oa the enor in flie wrapped phase map at JV^„ . 
Therefore, as the number of wavcjcixgttis used is 
increased the prooes3 reliability remams unchanged and 
the munber of ftinges that may be labelled Increases as 
gjvTO by cQuatiofl 6. 



Di^r^f LAvd Punafan i 



r 



mfWMm 



500 lOOD 

SDsnoloNomhar 



1 DiscFBte tcve) heterodyne j&mBtlons fer an optimam 4->. process. 

It should be noted that in the example of figure 1, the 
valna (f^jvy"* W ^ integer and therefore the discrete 
levels are equal hi phase values for each discrete level 
function. For any real positive number, 
algorithms ftinction correctiy- For ^axqple, 
//^^ =80.765 and using 4 projected fringe wavelengths 
the* discrete level fiitkCEions are given in figure 2. Tlie 
recursive imwrappmg process correctly unwraps 
successive beat finboigev tp obtata ftie J&Tooigc oxdsr, as seen 
in figure 3. 

Diseraia IbubI FuDctkn 1 



f — 

^ - ak ^..-..-ace* 


itt Sfto ._ t 1 e. «_3Cpo 


isao 


io 600 noo 




1900 



Fig. 2. PiwreUJ Jewel betcrodyiie^imcfloDsfer an optomn 4n^ process 
with atotal of j&lnscg. 




SatnpIffNUtnMr 

Fifi. 3. Ft'mge Or<ter MtoiJalol fbc «» ^mum 4-X prtcesa Tifiih a tawl 
or aQ.76S j^to^ 

For a typical four wavelength system a-^ » 2^r/85 and 
fee number of fringes that caa be conectjy idemti&ed witih 
6cr reliability is >1000, giving a dynamic range of 
>g5000. The following tdble shows values for Nj^ and 
dymwiic jMige for a range of X and cr^ . 



Phase 
Nbfse ' 
(rms, 
radians) 



2jr/l00 
2^r/20a 
2jr/A00 



Optimum 
3"^ Heterodyne 



Opiimmn 

4-X Heterodyne 

dynamic 
Range 



h 

13,095 
104J07 





163J00 

2,619,W0 
41,m800 



1,736 
13,889 

infill 

«88,8i59 

taWe I: Numbers of ftinges anfl dynamio rang© ottaincd with $1^.739% 
reUat>iniy for cptlmuxo 3 JL nnd 4-X hetenMlynB prooessing. 

il/terwctfiv^ M^Oi^hmfor Gefterating 
Optimum Wavelengths 

There altemativB ways to calculate the nimibers of 
fringe used at each wavelength that result in optimum 
oonfieuratlons In tiie discrete level functiQu^ For 

example with 3 wavelengdiB, //^^ , N^^ --^Nj-^ , 
Coixespondingly modi&cd expressions are 
then needed to form the heterodyne tractioxxs, equations 
4 and 5 may than be applied as before. It is expected ihat 
the noiae floor is worse in this case as it relies cn lihe 
difference of a greats nmnb^ of fringe sets. 
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3.L4 Modrficatiom to the Theoretical Description 
To aUow for es^>CTiniente]l eixox m obtaining exactly the 
number Qf jBilxges required across tiie measureniettt 
range, the wavelengflis would be selected close to their 
^eoretical values but in such a way as to define «n 
unambiguous measuremetit rasge sliglidy larger tiban tbfit 
jfBcpiibred. Therefore, in practical systems it would be 
cpmuontouse N^^.Nj^ -0,7, N^^-^ forattiree' 
wavelength system. 

4.0 Experimental Demonstration 

A wholle field, tdangulation based shape measurement 
systBio given as the ex^nplar for the theoretical 
analysis. This consists of a coherent fibre fiiage 
projector, wiiich produces a pattern of Young*a fringes 
across the test ofciject A CCD camera is then placed at an 
angle to the illumination dlrEctloti to view the fiinges, see 
fig. 4. Increasing the nvjmber of projected fiinge$ cr the 
angular sepajcatiop betweeii the CCD and the fiinge 
projector increases the sensitivity to the depth (z) of the 
test object. On© of the output fibres was wrapped aroimd 
a cylindrical PZT to allow ihe optical phase to be 
raodulated. An e^rfsting s«rvo control system allowed flie 
phase of ^e projected frmges to be stabilised by 
monitoring the reflections jBnom the distal ends of the 
output fibres the 4* ann of a directional coupler. The 
servo also enables accwaie 90"* phase steps to he 
obtained [11]. 




Flg.4. Experimavfi] setup 'fiv fiiase pip|ectiDn sfaape measurement 

A SBccMxd CCP camera was incorporated Into the fringe 
projector to measure the fibre separation, and hence the 
number of projected fringes. The fringes were sampled 
directly onto the second optical sensor via a polarising 
beamsplitter. A measurement of the phase distribution 
across the CCD gives the fibre reparation and hence the 
number of projected fringes [9], We Iwve d^aoaonstw^ted a 
resolution of <10 xun In the fibre separation measurement 
compared to a 50 mn requirement predicted by a 
shnulation for 3 projected wavelengths. 
To demonstrate the process, the flat side of an object was 
assessed initially as this gave easy verifioatiou of the 
fringe order calculation. Three wavelengths were used 
(from equation 3) with 100, 99 and 90 projected fringes. 



Fig. 5 rfiows the result of applyhig equations 4 and 5 to 
the three wrapped phase ro^9 obtained. Tlie central 
frtoge is awiomatlcaJly identified by ihe heterodyne 
process and is coloured black. Ihe fringe orders are 
clearly identified as a repaaitiing sdieme of 6 colours 
either side of the csgojtrai. ftioge. 




Fig. 5. EKvenmaM JSrtose map fer c^tiamm 3-X hetefodyoe 
processing (pscudccolour represcntaiion) 

The phase resolution obtained was estimated to be 1/SO* 
friQge pruxcipally Ihnited by the random speckle phase 
which contributes to the interfistencc phase. With 
tr^ ^27f/i0 and N^^=^109 equation 2 csaix be 
manipulated to determine the fringe order numberixig 
reliability e?cpBCted fron^ theory which evaluates to 
99.2%. A local heighbcmrhood check applied to the 
central flat area of the exijerimental data in fig- 3 shows 
the ttuaibep of pixels givimg the correct fringe order to be 
99.52%. Thereftteji a validation of the process and the 
underlying theory has been denaonstrated. 
The new process is equslly applicable to white light 
frioge pr^ijftctors based on spatial light modulators or on 
custom gratings produced on glass substrates. Compared 
to the reverse exponcntiaH temporal unwrapping scheme 
[8], to measure over -1€0 ftrnges requires 8 wrapped 
phase maps to be obtained Therefbre, the use of the new 
algorithm enables a rediiction data acquisition time 
and data space of >60H. The dyn^dic range obtained in. 
this measurement was 1:S000. FuitiiBnnorB, the process 
is applicable to any form of interferomeiric sensor where 
phase measurements are etvailable at defined sensitivities 
(wavelengths). 

5.0 Claims 

The following claims are loiade* 

1 , A method for absohrte fringe order calculation 
in multi-wavelength inierferometry. 

2. The wavelengths are selected in an optimum 
manner such that discrete level heterodyne 
functions contain equal numbers of levels in the 
interval -tt to -Hr. 
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A recursive unwrapping algoritiun has been 
fonnulated to determine ttie abiolute fHnge 
order ftoxn tibue discrete level functtOBS. 
The use of a phase noise measure (flxe 5te»dard 
deviation phase noise) provides ft tttcans to 
quantlftr reliabithy Of die fltige Older 

calculation. , j - 

A sot of mathematical rejanonships emhodymg 
claijns U 2 and 3 defining: the vavdengths to 
be wsed, a 9cbeme to fo«n the discrete level 
fiincdons, and a procedure to uttWtap these 
fhnctlons at each pixel independenily . 



3. 



4. 



5. 



<a) 



The number offiriDDgBsheterotb^ned to 6cr 
reliability is given by: 
2ff 



(b) 



6 A modification of claim 5 where the second 
wavelength equation <a» la modified by 

up » 30% in ord^ to guawntee foimtog ^ 
than 1 beat -fringe across the desired 
measurement ranger 

7. An alternative scheme fbr generating the 
optimum wavelengths and discrete level 
functions, which for a 3 wavdetjgth system are 

defined by N;^ . N^^ > V^"^ " 

8. Tbc process reliability may be defined to other 
values. For example, to obtain 8^ leliability, 
equatim (b) above becomes: 

^ 2g 

9. TTie process is appficable tp any foym of 
int^leromater, notjust to fHnge projection as 
exemplified herein* 



6.0 Essential Features of the luvention 

1 . A mechanism to produce a phase measurement at a 
number of discrete wavelengths. 

2. The wftY^engdi range required is dictated by 
equations (a) claim 5 or the equations in <slaim 7, 

3. A means by which the wavelengths (or tibie nranber 
of fi^inges across the measurement range) can be 
xneasured. 

4. A compnto* system to implement the data processing 
scheme and control of the data acquisition process. 



7.0 Naat-Coiifidcntial JJescription of flie 
iDventioTi 

This invention concerns the development of a roousi 
measmement process based on multiple wavelengftxs, 
phase measurement and heterodyne processmg. It is 
shown feat given a certain phase measurement ooise it is 
bossjble tp maximise the number of ftinge orders 
unambiguously identifiable and Hx^^^ ^ 
measurement dynamic range. Each fringe ^ ^ 
identified as an absolute immber, A robustness measure 
based on Gaufiaian statistics is also derived for the 
process. 
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